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We show that an LC parametric transducer can be effectively used to monitor an adiabatic 
evolution of the superconducting flux qubit. We propose a new scheme to measure the qubit's state, 
which is a quantum nondemolition measurement. The scheme can be easily extended to a three- 
qubit system, and allows the reading out of the qubits' states while the system remains in the ground 
state. An implementation of the adiabatic quantum algorithm MAXCUT for three superconducting 
flux qubits is discussed. 

PACS numbers: 85.25.Cp, 85.25.Dq, 03.67.Lx 



I. INTRODUCTION 

Ten years ago Peter Shor demonstrated theoretically-'^ 
that a quantum computer can factor large numbers much 
more effectively than a classical one. This discovery 
started an enormous effort to find a physical system 
which would be a suitable qubit, the building block of 
a quantum computer. Qubits are effectively two-level 
systems with controlled parameters. There are many 
systems in physics which can play the role of a qubit. 
One of them is a superconducting flux qubit which can 
be realized as a superconducting loop with low induc- 
tance Lq interrupted by three Josephson junctions. Its 
properties have already been analysed^'^ and experimen- 
tally verified.^ Superconducting qubits have several ad- 
vantages over qubits based on microscopic systems: they 
are scalable and can be accessed more easily and con- 
trolled individually. Moreover, aluminum technology, 
widely exploited for the preparation of conventional sili- 
con devices, can be used. 

Recently, several groups succeeded in demonstrating 
coherent macroscopic tunneling and Rabi oscillations 
in superconducting qubits. This can be considered 
as the first important step towards quantum computer 
realization.^"^ Most of them were time domain measure- 
ments, which are supposed to be important for quan- 
tum computing, since the much effort has been made in 
the direction of building a quantum computer based on 
a universal set of gates. However, in order to run an 
algorithm on such a universal quantum computer, quan- 
tum error corrections should be implemented. For a solid 
state qubit the error rate is only slightly below the thresh- 
old required for fault-tolerant computation. This places 
tremendous requirements on the hardware:^ the number 
of physical qubits should be larger than 10^ and tclcpor- 
tation between each two qubits should be possible. On 
the other hand, the new scheme of quantum computation 
based on adiabatic quantum evolution, which has been 
proposed by Farhi et al.,^ could solve tasks beyond the 
reach of present-day classical computers for a very moder- 



ate number of qubits 30). Very recently a scalable su- 
perconducting architecture for adiabatic quantum com- 
putation was proposed which requires nearest-neighbor 
coupling only.^° Moreover, Aharonov et al.^^ have shown 
that adiabatic quantum computation is equivalent to 
standard quantum computation. From an experimental 
point of view the adiabatic quantum algorithm MAX- 
CUT was demonstrated by an NMR technique on three 
qubit systems. In this paper, we propose a specific im- 
plementation for adiabatic quantum computing with a 
set of coupled superconducting flux qubits, which is pos- 
sible to realize with the present state of the art. We show 
that a parametric transducer can be effectively used to 
read out the results of the adiabatic evolution algorithm. 



II. PARAMETRIC TRANSDUCER AS A QND 
READOUT FOR ADIABATIC QUANTUM 
COMPUTATION 

Parametric transducers have been shown to be very 
sensitive instruments, that can overcome the standard 
quantum limit. The precision of the measurement of 
small changes of the dielectric susceptibility by a capac- 
ity transducer is of the order of 10~^°. In addition, a 
parametric transducer can work in a regime that satis- 
fies the criteria of quantum nondemolition (QND) mea- 
surements. Usually, an electromagnetic auto-oscillator is 
used as a key element of a parametric transducer since 
the frequency can be measured with a very high accuracy. 
The scheme of a parametric transducer is shown in Fig. I, 
it contain a high quality LC resonator connected to an 
amplifier. The resonant frequency of the LC circuit de- 
pends on both the inductance L and the capacitance C 
by the relation ujr = 1/VLC. In our experiments, typi- 
cally LUr /2t: ~ 30 MHz. This satisfies uJr <^ Wq, where Uq 
is the transition frequency between the ground and first 
excited energy level of the qubit. Thus, the magnetic 
susceptibility of the qubit placed in a resonator can be 
measured from the shift of the resonance frequency. It 
can be easily shown^'^ that the tangent of the phase shift 
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FIG. 1: Scheme of a parametric transducer inductively cou- 
pled to a superconducting flux qubit. The rf voltage across 
the tank is amplified by a cooled HEMT amplifier thermally 
linked to a f K pot.^* After room temperature amplification 
the signal is detected by an rf lock-in voltmeter. Both the 
amplitude and phase of the r f voltage are measured as a 
function of the external magnetic fiux applied to the qubit 
produced by the currents Idc and lb through a coil and a 
wire, respectively. 



9 between the voltage across the tank and driving current 
is proportional to the real part of ac susceptibility x' 
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-k'Qx' 
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where < A: < 1 is the coupling coefficient between the 
resonator and sample. The ideas behind a parametric 
transducer were also used in the design of an r/-SQUID 
by Silver and Zimmermann.-'^^ It was shown theoreti- 
cally that an r/-SQUID can achieve the quantum limit. 
Therefore, the parametric transducer is a suitable read- 
out device for superconducting flux qubits. 

The magnetic susceptibility of the superconducting 
flux qubit is"'^^ 



A2 



(A2 -(-£2)3/: 



■ tanh 



T 
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where A is the tunneling amplitude, Lq is the induc- 
tance of the flux qubit, Iq is the persistent current in the 
qubit, T is the temperature, and e = ^olqf is the bias 
of the qubit, where / is the deviation from degeneracy 
defined in terms of internal magnetic flux in the qubit as 
/ = ^i/^o - 0.5. By using Eqs. (1),(2) the persistent 
current and the tunneling amplitude can be determined 
experimentally by measuring the resonator phase as a 
function of the external magnetic flux $e-^^ The func- 
tion x'if) (Eq. 2) has a simple form, and it is easily seen 
that x'if) exhibits a peak at the degeneracy point f — 0. 
If the temperature T ^ A, the explicit equations for the 
persistent current and the tunneling amplitude can be 
readily derived 
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where x'a ^^'^ /fwhm £^re the peak amplitude and the 
full width at half maximum (FWHM), respectively. 




FIG. 2: The phase shift 6 between the bias current 7^/ and the 
rf voltage of the parametric transducer inductively coupled 
to the superconducting flux qubit as a function of the internal 
magnetic fiux in the qubit. The curve with hysteretic behavior 
(black curve) corresponds to the 'qubits' with a large ratio g — 
Ej/Ec ~ 10^ (classical regime). The straight line (vertically 
shifted for clarity) and the non-hysteretic line correspond to 
qubits with g ~ 60, and a — 0.9 and a = 0.8, respectively. 



Here we would like to point out that the measure- 
ment by means of a parametric transducer is a quan- 
tum nondemolition measurement, because the qubit is 
staying in its ground state the entire time of the mea- 
surement, as the resonant frequency of the resonator cUr 
is much lower than the transition frequency LUq. The 
output signal of the parametric transducer contains in- 
formation about the amplitude of the persistent current, 
but holds no information about the phase of the rapidly 
oscillating persistent current. A parametric transducer 
cannot even distinguish whether the current flows clock- 
wise or counterclockwise. This can be directly seen for 
'qubits' in the classical regime (see the hysteretic curve 
in Fig. 2). Exactly at the degeneracy point f — 0, the 
two branches of the hysteretic curves corresponding to 
current flowing clockwise and counterclockwise cross, i.e. 
the transducer gives the same signal. The reason for this 
is that the operator probed by the inductive transducer 
is ax as we show below. In this sense, such a readout is 
complemental to the SQUID readout which measures tr^ 
(ax and cr^ are Pauli matrices) to make the analogy with 
Stern-Gerlach apparatus complete. More formally, the 
Hamiltonian of a qubit-resonator system at the degener- 
acy point / = can be written in the form2'^ 

H = Hr + Hg + Hint = hujrb^b + Aax + j{b^ + h)az (5) 

where 5^,6 are creation and annihilation operators, re- 
spectively, of the photon field in the resonator, 7 — 
k^fiLiJrLqlq is the coupling energy between the resonator 
and qubit, and ax and az are Pauli matrices of the nat- 
ural basis of the qubit (i.e. the two eigenstates of opera- 
tor az correspond to the currents fiowing clockwise and 
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counterclockwise). After unitary transformation 

the Hamiltonian (5) takes the form 

UiHUl = hwrb^b + Act; + 7(6^^ + ba'+) (7) 

where 
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are spin-flip operators and cr^ is the Pauli matrix in the 
eigenbasis of the qubit at the degeneracy point. Following 
the approach in Ref. 21, after a second unitary transfor- 
mation 

C^2=exp(^(6a;-6V)), (9) 

and by expanding to second order in 7/ A, the trans- 
formed Hamiltonian H' = U2U-iHUlU2 is 



The readout procedure could be as follows; let us sup- 
pose that the qubit is in the state |1) (i.e. / > 0, see 
Fig. 3). If the internal magnetic flux in the qubit is 
changing towards zero, then the qubit is moving through 
its degeneracy point (/ = 0) where two classical energy 
levels cross (dashed lines in Fig. 3). At this point the 
qubit is in the superposition of the states |0) and |1) 
where the magnetic susceptibility of the qubit changes 
rapidly. Thus, the inductive transducer gives a consid- 
erable signal. On the other hand, if the qubit is in the 
state |0) (i.e. / < 0), one should increase the external 
magnetic flux to move the qubit to the degeneracy point. 
If we do not know the state of the qubit at / = /o we 
can sweep the external magnetic flux in order to change 
the internal magnetic flux in the qubit around this point 
and from the response of the parametric transducer we 
can determine whether the qubit was in state |0) or |1) 
(signal is observed for / > /o or / < /o, respectively). 
In the next section we will show numerically that the 
qubits can be readout one after another while staying all 
the time in the ground state of the system. 
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where Wq — Lql'^/2 is the magnetic energy of the qubit. 
Experimentally, a shift of the resonant frequency of the 
resonator, which depends on the qubit state, is mea- 
sured. This shift is determined by the first term of Eq. 10, 
i.e. the measured observable is cr^, and one can readily 
find that the sufficient condition for QND measurement 
[(7^, H'] = is satisfied. -'^■^ Provided that the coupling be- 
tween resonator and qubit is small a'^ corresponds to 
in the original basis. This means that the resonator mea- 
sures the observable at the degeneracy point in con- 
trast with the SQUID which measures cr^. Let us point 
out additional difference between resonator and SQUID 
measurement. The SQUID measurement makes a pro- 
jection of the spin into the 2;-axis, i.e. at the degeneracy 
point the qubit is localized in one of the classical states 
after the measurement, and the SQUID gives a signal 
corresponding to this state. This measurement is non- 
QND since the SQUID is coupled directly to the oscillat- 
ing variable. On the other hand, the resonator gives no 
signal if the qubit is in the eigenstate of the operator Uz, 
i.e. the resonator does not perform a measurement, and 
therefore, does not disturb the qubit. Such a readout 
method has a clear advantage in the case of adiabatic 
quantum computing. The qubit remains in its ground 
state also after the measurement, i.e. the measurement 
of one qubit does not spoil the result of the adiabatic evo- 
lution. However, it should be noted that this statement 
is valid only if the amplitude of the circulating current 
in the resonator is small enough to avoid Landau-Zencr 
transitions. Nevertheless, as we have shown theoreti- 
cally in Ref. 23 and experimentally demonstrated in sec- 
tion IIIB, the noise temperature of the cooled amplifier^^ 
enables one to fulfill this condition. 



III. ADIABATIC EVOLUTION 
A. Theory 

The idea of quantum computation by adiabatic evolu- 
tion is very simple but, surprisingly, was discovered only 
recently.^'^° It is based on the fact that, in practice, it is 
very difficult to find a ground state of certain Hamilto- 
nians. Such a task belongs to the set of non-polynomial 
time (NP) problems. On the other hand, some Hamilto- 
nians have a trivial ground state which is easy to find. 
Let us assume that the Hamiltonian of N qubits H{p) 
can be externally controlled by the parameter p and that 
its ground state is separated from the first excited state 
by the energy gap g{p) = Ei{p) - Eo{p) {Eq, Ei are the 
two lowest eigenvalues of the Hamiltonian H{p)). Pro- 
vided that the ground state of the initial Hamiltonian 
Hj = H{p = 0) can be easily found, we can construct it 
and then change the parameter p slowly from p = to 
p — 1. If we do it sufficiently slowly, i.e. in a time t ^ 
^Smax/gmin where the Emax ~ ma,xEi{p) - mmEo{p) 
and gmin = ^^^9{p), the ground state of Hj is evolved 
to the state which is with high probability the ground 
state of Hp = H(p = 1). Thus, we have prepared the 
system of the qubits in the ground state of Hamiltonian 
Hp and they can be read-out. As a matter of fact the sys- 
tem is in the ground state of the Hamiltonian H{p) dur- 
ing the whole adiabatic evolution, i.e. the system is im- 
mune against dephasing and relaxation. Here we should 
emphasize that the adiabatic evolution of the Hamilto- 
nian is crucial in speeding up considerably the finding 
of the ground state of the Hamiltonian Hp. One could 
suggest that it is enough to wait a while and the sys- 
tem would relax itself into the ground state. However, 
a Hamiltonian which encodes a NP problem exhibits a 
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lot of local minima and the physical system needs an 
exponentially long time (as a function of the number of 
qubits) to find its global minimum. As an example, one 
can consider Ising model of N antiferromagnetically cou- 
pled magnetic moments. It is well known that such a 
system can be highly frustrated. The task of finding the 
minimum of the Ising Hamiltonian is equivalent to the 
optimization MAXCUT problem which belongs to a NP- 
complete problem.^'* Thus, it seems that NP problems 
cannot be solved in polynomial time on either digital or 
analog classical computers. Theoretically it was shown^^ 
that an adiabatic quantum algorithm can find the global 
minimum of some functions in polynomial time whereas 
a classical simulated annealing algorithm requires expo- 
nential time. The crucial condition for adiabatic quan- 
tum evolution is the existence of an energy gap between 
the ground and upper levels. This is the key difference 
between classical and quantum systems, thereby enabling 
an enormous speed up of adiabatic quantum algorithms 
over classical ones. The size of the energy gap limits the 
speed of adiabatic quantum evolution as we will show 
experimentally in the next section. 

Adiabatic evolution can be demonstrated on a single 
qubit. Following the original paper by Farhi et al.,^ we 
start from the initial Hamiltonian at i = 



Hr 



(11) 



Then we adiabatically evolve from Hj to the problem 
Hamiltonian Hp in time r 



Hf 



(12) 



This scheme can be implemented for a superconducting 
flux qubit. Near the degeneracy point / = 0, the qubit 
can be described by the Hamiltonian 



H{t) = e{t)a^ + Aoi 



(13) 



At a bias e = 0, the two lowest levels of the qubit anti- 
cross (Fig. 3), with a gap of 2A. By increasing e slowly 
enough, the qubit will adiabatically transform from the 
superposition state (|0) + |l))/\/2 to but remains in 
the ground state. For \s{t)\ ^ A, A diminishes and the 
Hamiltonian takes the form 



H{t) = e(T)oi 



(14) 



However, if the bias changes in time e{t) — Xt, the 
qubit can 'jump' from the ground state \g) to the ex- 
cited state |e) with probability Plz — exp(— 7rA^/?iA). 
This process, known as a Landau-Zener transition, ■^^ 
would violate adiabatic evolution and, therefore, should 
be avoided. This puts constraints on the characteris- 
tic time T of the adiabatic evolution which can globally 
be estimated as: r » UEj/A^. Consequently, r can 
be considerably shorter if we take into account that a 
Landau-Zener transition takes place only in the A vicin- 
ity of the anti-crossing point. Thus, e(t) can be changed 
quickly except in the region close to anti-crossing point. 




FIG. 3: Quantum energy levels of the qubit as a function 
of normalized internal magnetic flux/ = $i/$o — 0.5. For / 
much less or greater than zero, the qubit is in the state |0) or 
1 1), respectively. The dashed lines correspond to the classical 
potential minima. 



For such a local adiabatic evolution the requirement for 
T reads r 3> h/ A. Note that only this condition leads to 
a quadratic speed-up of the adiabatic evolution version 
of Grover's algorithm. A measurement by a parametric 
transducer provides the unique possibility of controlling 
the speed of an adiabatic evolution. The smaller the en- 
ergy gap is the larger is the signal from the transducer 
(see Eq.1,2). This signal can be used as feedback for e(t) 
sweeping so that the condition for adiabatic evolution 
can be satisfied locally for an unknown ground state of 
the system. 

The tunnel splitting 2A is very sensitive to the Joseph- 
son and Coulomb energy of the junctions. It can be finely 
tuned by reducing the size of one junction in the super- 
conducting loop, while leaving the two others unchanged. 
If the ratio between the area of the small and large junc- 
tion is a (a < 1), A can be roughly estimated^ ■^'^ 
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where g — Ej/Ec- By changing the parameters a and 
5, one obtains a crossover from the classical, through the 
Landau-Zener, to the adiabatic regime. 



B. Experiment 

In order to demonstrate the crossover from the classical 
to quantum regime we have prepared three qubits with 
different parameters a and g. The qubits were placed 
inside pancake niobium coil made by using electron- 
beam lithography on oxidized Si substrates. The typi- 
cal linewidth and the distance between 20-30 coil wind- 
ings are 1-2 /xm. The coils self-inductances are L — 
50 — 140 nH. For all experiments reported here we use an 
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external capacitance C = 470 pF, therefore the tank res- 
onance frequency is 19.6 — 32.8 MHz with quahty factors 
Q=700-1700. The 3JJ qubit structure was fabricated out 
of Al in the middle of the coil by a conventional shadow 
evaporation technique. The critical current was deter- 
mined, by measuring an rf- SQUID prepared on the same 
chip, as Ic — 250 — 400 nA. The qubit's loop area was 90 
/im^, with Lq = 40 pH. The typical coupling coefficient 
between the coil and qubit is 1 — 2 x 10"^. In Fig. 2 the 
typical response of the inductive transducer is shown for 
three values of the parameters a and g, which correspond 
to three different regimes: classical, Landau-Zener and 
adiabatic. In classical regime the signal from paramet- 
ric transducer is proportionate to the first derivative of 
the Josephson current with respect to internal magnetic 
flux.'^^ Close to degeneracy point there are two classi- 
cal states corresponding to the currents flowing clockwise 
and counterclockwise (hysteretic behavior). For g — 
and a — 0.9 the qubit is in an intermediate regime where 
both tunneling between two classical states and Landau- 
Zener transitions are not negligible. There is still no vis- 
ible dip in the phase characteristic but the losses caused 
by Landau-Zener transitions decrease the quality factor 
of the resonant circuit and, consequently, the amplitude 
of the rf voltage. By keeping g constant, but decreasing 
the size of the third junction from a — 0.9 to a = 0.8, the 
tunnel splitting 2A increases and Landau-Zener transi- 
tions are suppressed. As a result, a shift of the resonance 
frequency of the parametric transducer leads to huge dips 
in the vs f curves (adiabatic regime). Nevertheless, if 
the voltage amplitude across the parametric transducer 
is increased high enough, the Landau-Zener transitions 
suppress the dip again. Under this condition, a discrep- 
ancy between experimental and theoretical curves, calcu- 
lated within the adiabatic approach is observed (Fig. 4). 
Thus, we have observed the crossover from the classical, 
through the Landau-Zener, to the adiabatic regime of a 
superconducting flux qubit by decreasing the size of the 
Josephson junctions. Our experimental results show that 
the idea of adiabatic quantum computing can be demon- 
strated on a system of superconducting flux qubits. A 
reasonable and primarily feasible design is shown in the 
next section. 
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FIG. 4: The phase shift 6 between the bias current /,./ and 
the rf voltage of the parametric transducer as a function of 
the normalized internal magnetic flux for small Vrf ~ 
(lower curves) and large Vrf ~ 5/iV (upper curves) rf volt- 
ages. The resonant frequency of the parametric transducer 
was 32 MHz. The discrepancy between experimental (solid 
line) and theoretical (dotted line) curves for the large ampli- 
tude rf voltage is caused by Landau-Zener transitions. 



can be encoded into a Hamiltonian H of N inductively 
coupled superconducting qubits 



H 



N 



N 

1 = 1 



N 

i<j 



i.j^Zji^z.j 



(17) 



where and (Tz are Pauli matrices, ei{fi) is the energy 
bias of the z-th qubit, and is the coupling energy 
between the i-th and j-th qubit. The eigenvector |s), 
corresponding to the ground state of the Hamiltonian 
H, is the solution of the payoff function P{\s)) if (a) 
^ Jij Vi,j, and (b) Si = —Wi/2, and J^^ — Wi^j/2. 
For superconducting qubits, the initial Hamiltonian Hj 
can be easily constructed by taking into account that 
J,,j = and A, = if /i = -0.5, i.e. 



Hj 



N 

E^ 

i=i 



.(-0.5)(T, 



(18) 



IV. IMPLEMENTATION OF THE MAXCUT 
PROBLEM FOR A SET OF INDUCTIVELY 
COUPLED SUPERCONDUCTING QUBITS 

The MAXCUT problem is a part of the NP-complete 
problems. Mathematically, in order to solve the MAX- 
CUT problem, one should find the maximum of the payoff 
function^^ 



(16) 



where Wij , Wi are the parameters of the problem and 
Si = 0, 1 are components of the vector \s). The problem 



The ground state of Hi is trivial, |0). By changing the 
bias of individual qubits adiabatically to Si — —Wi/2, 
the Hi is transformed to H . (The coefficients Wij are 
set by design and they are determined by coupling en- 
ergies between qubits.) H encodes the payoff function 
P{\s)) completely if A; = 0. Unfortunately, we can- 
not switch off the tunnel splitting A^ in superconducting 
qubits, but it is not absolutely necessary if Ji^j ^ A^. 
Nevertheless, we will show that by making use of a para- 
metric transducer, one can obtain the answer even if 
Jij ^ Ai. Moreover, the qubit states can be readout 
while staying in the ground state of the system. 

The most simple but still reasonable example of the 
adiabatic quantum optimization algorithm MAXCUT 
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FIG. 5: Three-Qubit design for MAXCUT problem, (a) 
Three superconducting flux qubits are placed in supercon- 
ducting coil. The qubits can be biased independently by 
dc bias wires 61,62,63. (b) The qubits are coupled through 
common Josephson junction marked by circles. The coupling 
energy'^" Jij = {Mij +$o/2TvIc)IqiIqj, where Ic is Josephson 
critical current of the common junction. 



can be implemented by three coupled superconduct- 
ing flux qubits {N = 3). The coupling between the 
qubits can be reahzed by means of a common Joseph- 
son junction'^*' shared between two qubits (see Fig. 5). 
This enable to increase the coupling energy over pure 
magnetic one. The couphng energy « 0.3 K has been 
measured recently.^^ Thus, for the present design we have 
chosen the interaction energies between the qubits to be 
Ji^2 = J2.3 — J1.3 — 0.3 K, the persistent currents are 
Ipi — Ip3 — 350 nA, Ip2 = 420 nA and tunneling ma- 
trix elements are Ai = A2 = A3 = 96 mK. By choos- 
ing appropriate values for it is possible to realize the 
situation that the system exhibits both a local and a 
global minimum. We have chosen the following param- 
eters ei(0.006) = 0.315 K, £2(0.004) = 0.252 K, and 
£3(0.01) — 0.525 K. The energy of the ground state for 
various vectors \s) is shown in Table L In the state |101) 
the system is in the global minimum. Note that for |110) 
the system exhibits a local minimum, that is, there is no 
way to decrease the energy of the system by flipping the 
persistent current in one qubit only. Thus, the system 
can stay in the state |110) for an exponentially long time 
at low temperatures. In our design the lowest 'energy' 
barrier which the system sees from the local minimum is 
higher than 0.5 K. This could lead to a wrong answer, 
unless the Hamiltonian transform is carried out adiabat- 
ically. 

The qubits' state can be readout by an inductive trans- 
ducer as was described above. The internal magnetic 



TABLE I: Energy of the system for various vectors. Ji,2 = 
J2.3 = Ji,3 = 0.3 K. El = 0.315 K, £2 = 0.252 K, and £3 = 
0.525 K. 



flux of the individual qubits can be changed by a current 
through the wires placed nearby each of them. In such a 
configuration, all three qubits can be readout by making 
use of one transducer only. Nevertheless, the idea should 
be checked since qubits interact and A is nonzero. The 
three qubit Hamiltonian can be solved numerically. In 
the following section we simulate the readout of a para- 
metric transducer inductively coupled to three supercon- 
ducting flux qubits. 



A. Numerical simulation 

The Hamiltonian (17) was solved numerically and the 
energy levels of the Hamiltonian (17) as a function of fi 
are shown in Fig. 6. We have used the same parame- 
ters as those used in our design. We have also calcu- 
lated the response of the parametric transducer using the 
formula'^^^'^^ 



tan 9 



-2Qy— ^ 



(19) 



where E^, — Eq is the distance between ground and upper 
energy levels, and 

Rau = |^^fc,/L~/,,(0|a«|z^>^ 

X |^E^^V^^«(H^P^|0)j (20) 

are the real matrix elements. Here fc; is the coupling co- 
efhcient between the i-th qubit and resonator, and Lqi 
is the qubit's inductance. For Q = 1000, L = 81 nH, 
Lqi — 40 pH, and ki = 0.036 the results are shown in 
Fig. 7. From these figures it is apparent that the qubits' 
states can be determined by a parametric transducer. We 
have also tried to find the threshold for A^ below which 
the state of the qubit cannot be distinguished. As a cri- 
terion the existence of the distinguishable dips on the 
experimental curves can be chosen. From Fig. 8 one can 
see that the positions of the dips do not change as A in- 
creases and they can be distinguished for relatively large 
values of A. Thus, the parametric transducer readout 
delivers the right solution of the problem. 



V. CONCLUSIONS 

Experimentally, we have demonstrated the principle of 
adiabatic quantum evolution in a single qubit. Theo- 
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FIG. 6: First three energy levels of the three qubit system 
during readout. Readout of the qubit starts at point ff — 
0.006, /f = 0.004, and /f = 0.01 then its bias is changed 
adiabaticaly and separately through qubit 1 (a), 2 (b), and 3 
(c) while keeping it fixed in the others. At the points with a 
large curvature of the ground level the parametric transducer 
gives a considerable response (see Fig. 7). If this point is on 
the left (right) side of the point corresponding to the problem 
Hamiltonian (marked by green vertical line) the qubit is (or 
better to say would be if A = 0) in the state |0) (|1)). 



retically, we have shown that three inductively coupled 
superconducting flux qubits placed in a superconducting 
coil can be used to demonstrate the adiabatic quantum 
algorithm MAXCUT which belongs to the set of NP- 
complete problems. A three qubit design has been pro- 
posed and simulated numerically. 

Note added in proof: Recently, Lupascu et al. [A. 




-1.0 -0.5 0.0 0.5 1.0 

FIG. 7: The phase shift between the voltage and bias current 
of the parametric transducer with respect to A/ = / — f^. 
Readout of the qubit starts at point /f — 0.006, /f = 0.004, 
and /g = 0.01 then its internal magnetic flux is swept adiabat- 
icaly around this point. The red (solid), black (dotted) and 
blue (dashed) lines correspond to bias flux change in qubit 1,2 
and 3, respectively. From the position of the dips we find that 
the state |101) corresponds to the global minimum (compare 
with Table I). 



0.0 F 




-1.0 -0.5 0.0 0.5 1.0 

FIG. 8: The phase shift between voltage and bias current 
of the inductive transducer with respect to A/ = f — 
for various values of Ai (Ai is taken to be the same for all 
qubits). Readout of the qubit starts at point /f = 0.006, /f = 
0.004, and /f — 0.01 then the internal magnetic fiux is swept 
adiabaticaly around this point. The red solid line and black 
dotted line correspond to the qubits 1 and 2, respectively. 
From the upper to lower curve (at A/ = 0.01) Ai takes the 
values 0.048, 0.096, 0.144, 0.192, and 0.240 K. 



Lupascu, C. J. M. Verwijs, R. N. Schouten, C. J. P. M. 
Harmans, and J. E. Mooij, Phys. Rev. Lett. 93, 177006 
(2004)] proposed a similar readout method which enables 
to measure the observable tr^ in a nondestructive way. 
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